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ABSTRACT 


The first numerical study of axial (toroidal) pulsation modes of a slowly rotating relativistic 
star is presented. The calculation includes terms of first order in e = Sly/R 3 /M « 1 (R is 
the radius, M is the mass and SI is the rotation frequency of the star), and accounts for effects 
due to the coriolis force. Effects due to the centrifugal flattening of the star enter at order e 2 
and are not included in the analysis. It is shown that increased rotation tends to decrease the 
damping times for prograde modes, while retrograde become longer lived. Specifically, we show 
that rotation affects the axial gravitational-wave w-modes in this way. We also present the 
first relativistic calculation of the so-called r-modes (analogous to Rossby waves in the Earth’s 
oceans). These have frequencies of the same order of magnitude as the rotation frequency of 
the star. The presented results indicate that the r-modes are unstable due to the emission of 
gravitational radiation for all rotating perfect fluid stars. This is interesting since the previously 
considered gravitational-wave instability associated with (for example) the /-mode of the star 
sets in at a critical rotation rate. Because they are unstable also for the slowest rotating stars 
the r-modes may well be of considerable astrophysical importance. 
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1. Introduction 

This paper is the first in a series of investigations of perturbations of slowly rotating 
relativistic stars. The original intention was to present an exhaustive description of the problem, 
together with a detailed discussion of the effect that rotation has on the various pulsation 
modes of a relativistic star. However, a preliminary investigation into the axial (also referred to 
as toroidal, or odd-parity) mode-problem unveiled an interesting and rather surprising result: 
There exist pulsation modes of a rotating relativistic star that are unstable due to the emission 
of gravitational radiation at all rates of rotation. That is, these modes — the relativistic 
analogue of the Newtonian r-modes, cf. Papaloizou and Pringle (1978) and Saio (1982) — are 
unstable even for very slowly rotating perfect fluid stars. The instability of the r-modes is thus 
different to previously considered mode-instabilities (e.g. for the /-mode of the star) that set in 
at a certain rate of rotation (Comins 1979; Friedman 1983; Managan 1985; Imamura, Durisen 
and Friedman 1985; Lindblom 1986; Lindblom and Mendell 1995; Stergioulas and Friedman 
1997). But the fact that the r-modes are unstable can still be understood in terms of the 
mechanism that was first discussed by Chandrasekhar (1970) and Friedman and Schutz (1978). 

This paper describes the calculation that led to the discovery of the unstable r-modes 
without many of the technical details. These will be reported elsewhere. 


2. Slowly rotating relativistic stars 

This paper concerns pulsation modes of a slowly rotating relativistic star. The modes follows 
from a study of linear perturbations of the stellar fluid and the associated spacetime metric, for 
a star that rotates sufficiently slowly that the effect of rotation can be considered as a small 
perturbation of a non-rotating configuration. The analysis is based on a double perturbation 
expansion. To lowest order the star is static and non-rotating, and the perturbations are such 
that \h fJ/U \ << 1 , where g = ^background _|_ and a i so £ = Q^/R3/M « 1. Here R is the 
radius of the star, M is the mass, and f 1 is the rotation frequency of the star according to a 
distant observer (the star rotates uniformly so fi is a constant). Geometrized units c = G = 1 
are used throughout the paper. 

Before presenting the details of the problem it is meaningful to discuss whether a “slow- 
rotation expansion” is likely to provide results of physical relevance. It is rather straightforward 
to show that this is the case. The rotation of a star is absolutely limited by the Kepler frequency, 
at which mass shedding at the stellar equator makes the star unstable. The following empirical 
formula has proved to be a reasonable approximation of this limiting frequency 



where C « 7.8 x 10 3 s -1 (Haensel and Zdunik 1989; Friedman, Ipser and Parker 1989; Cook, 
Shapiro and Teukolsky 1994). M and R refer to the mass and the radius of the corresponding 
non-rotating star. By rewriting the definition of the expansion parameter e in the spirit of this 
relation we get 
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where D = 8.6 x 10 -5 s. That is, the limiting Kepler frequency corresponds to a value of e 
that is slightly smaller than 0.7. This is the largest value of e that should ever be relevant. A 
calculation including only terms of order e will probably not provide an accurate description of 
extremely rapidly rotating stars, but it is interesting to note that an approach that is consistent 
to 0(e 2 ) may prove useful also for such cases. It has, in fact, been shown that a second-order 
calculation (e 2 ) of the upper mass limit (and other properties) for rotating neutron stars is 
accurate to within a few percent (Weber and Glendenning 1992). 

As a second example, consider a millisecond pulsar with a period of 2 ms (cf. PSR 1937+21 
that has a period of 1.6 ms). This corresponds to fl = 3100s -1 and if we for simplicity assume 
that R = 10 km and M = Mq the corresponding value of e is roughly 0.3. This value should 
be within reach of a calculation to first order in the rotation expansion. Consequently, such 
a calculation has the potential to provide considerable information about the fastest spinning 
pulsars that have been observed. 

There are practical reasons why the perturbation approach may be preferred to a fully 
nonlinear calculation. There are by now several reliable methods for numerically solving the 
Einstein equations to construct fully relativistic rotating stellar models (see, for example, 
Friedman et al. 1986; Stergioulas and Friedman 1995), but to calculate the neutral modes 
(that signal the onset of gravitational-wave instability) for these models is still a difficult task 
(Stergioulas and Friedman 1997). And it will be considerably harder to calculate the actual 
pulsation modes of a rapidly rotating, relativistic star. A full consideration of the pulsation 
properties of rapidly rotating stars will require new computational techniques and considerable 
computer power. In contrast, the slow-rotation approach should not require any conceptually 
new steps. The relevant perturbation equations can be reduced to wave equations, and the 
main difference from the non-rotating case is that one must consider several coupled equations. 
That the equations take a standard form is an enormous advantage when one is interested 
in extracting the complex frequencies of the pulsation modes. Moreover, the slow-rotation 
calculation serves as a useful benchmark test for future work on rapidly rotating stellar models. 


3. Axial perturbations of a slowly rotating star 


As was shown by, for example, Hartle (1967), the centrifugal force affects the shape of a 
rotating relativistic star through terms of order e 2 . Hence, one can assume that a slowly rotating 
star remains spherical as a first approximation. The corresponding metric can be written 

ds 2 = —e"dt 2 + e x dr 2 + r 2 (dd 2 — 2wsin 2 ddtdip + sin 2 9d<p 2 ) (1) 


For a given equation of state p = p{p), where p is the energy density, the pressure p and the 
two metric coefficients v and A are determined by the standard TOV equations [I use standard 
notation and the relevant equations are listed by, for example, Kojima (1992)]. 

When terms of order e are included in the slow-rotation expansion one must account for the 
“dragging of inertial frames”, that is represented by u> in ( 0 )- After defining u> = Q — oj we have 
(Hartle 1967) 


d 2 oj 

dr 2 


47 r(p + p)e x r 


4 dO 
r dr 


167r(p + p)e x ui = 0 . 


( 2 ) 
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The desired solution to this equation is well behaved at the centre of the star; 


1 w 0 


1 + ^r(Po + Po)r 2 + 0(r 4 ) 
5 


as r 


0 


( 3 ) 


where po and po are the central values of the pressure and the density. In the vacuum outside 
the star we have 

" ( 4 ) 


U) = 

where J is the total angular momentum of the star. The constants wo and J are fixed by the 
requirement that w, and its derivative, be continuous across r = R. 

The equations that describe perturbations of a slowly rotating star in general relativity have 
previously been considered by Chandrasekhar and Ferrari (1991a) and Kojima (1992). I will use 
the notation of Kojima (1992), and refer the interested reader to his paper for a listing of the 
original perturbation equations. 

The perturbations of a non-rotating star can be divided into two classes: axial and polar 
(often referred to as odd- and even parity perturbations, respectively). The equations that 
describe these two classes decouple and can be considered separate from each other. Specifically, 
one can show that axial perturbations are governed by a single wave equation (in Regge-Wheeler 
gauge) 

d 2 X, 


Vf(X lm ) = 


L Zm 


drl 


+ (a 2 - V l x )X lm = 0 . 


( 5 ) 


Once the function Xi m is known one can infer the metric perturbations hte, ht v , h r g and h rt p 
as well as the fluid displacement £‘0 The angular dependence of all perturbed variables is 
expressed in terms of the spherical harmonics Y[ m [for details, see Kojima (1992)]. In (|j) it is 
assumed that the time dependence of the perturbation is e~ lat . The assumption of a harmonic 
time-dependence is adequate as long as the focus is on the spectral properties of the star. The 
tortoise coordinate r*, which is defined by 


_d_ 

dr * 


= e 0-A)/2 


dr 


is familiar from studies of perturbed black holes, and the effective potential is 


V l x (r)=e" 


l (l 


1) 6M , . 

72 - -^T+Mp-p) 


( 6 ) 


( 7 ) 


This expression is valid both inside the star and in the exterior vacuum [M(r) is the mass inside 
radius r]. 

Let us now assume that the star is slowly rotating. To incorporate the first order effects of 
rotation on the perturbations we consider 


Xlm=X^ 


X, 


(1) 

Im 


( 8 ) 


where xj^ is a solution to the axial equation for a nonrotating star and xj^j corresponds to 
terms of order e. This results in the coupled equations 


= 0 , 


and 


V? (I ; W) = At + Vtlm + V?-i m ■ 


( 9 ) 


( 10 ) 
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In this schematic description, Af rn represents the coupling to the zero order axial equation, and 
V^ ±lrn couple the first order axial perturbations to the zeroth order polar ones. 

The explicit form of the coupling to the lowest order axial solution is (Kojima 1992) 


_4; y - 


2m 


1(1 + l)a 

z u duj 


{ [l(l + l)tr 2 (fl — w) + [1(1 + 1) — 2]87r(p + p)e 1 'u> 


— -^-—{[/(/ +1) — 2] (2r — 5M — 47rpr 3 ) — 2re A ) 

r z dr 


A, 


(0) 

lm 


i(i + i )e (— A )/ 2 —^ 

dr dr 


(0) 


( 11 ) 


In the exterior vacuum this expression reduces to 


A* 
^lm 


2 mj 
1(1 +l)c 


61(1 + \)e v dX\ 
r 4 dr * 


2Z(Z + l)cr 2 _ [Z(i + l)-3](2r-5M)-M 


X, 


(o) 

lm 


( 0 ) 


( 12 ) 


The explicit form for V(± lrn was given by Kojima (1992). 

The apparent coupling between the axial and polar perturbations in (10) potentially 
complicates the definition of a stellar pulsation mode, but these complications enter only at 
order e 2 and higher. As in the non-rotating case, the asymptotic behaviour of a general solution 
is 


Xlm ~ XoutC + XmC 


(13) 


for a general frequency cr. Here x ou t = x^ t + x ^ t , and similar for Xi n . In analogy with the 
non-rotating case a quasinormal mode of the star should correspond to a frequency a n such 
that x- ln (a n ) = 0. In principle, one can find such solutions to the axial slow-rotation equations 
for any magnitude of 'Pf^_ lrn (Chandrasekhar and Ferrari 1991a). But in this paper the primary 
interest concerns modes that limit to the axial modes of a nonrotating star as e —> 0. For such 
modes the zero order (nonrotating) polar perturbations should vanish, and hence one should use 
V^ +lrn = ViL lm = 0 in ([To]). The coupling between the two classes of perturbations will affect 
the mode-frequencies only at order e 2 and higher (Kojima 1993a). In principle, an axial mode 
generates a polar perturbation of order e that feeds back and affects the axial mode-spectrum at 
second order. Since this paper is restricted to the first order axial mode-spectrum the effect of 
the coupling to the polar perturbations will not be included. An axial pulsation mode is defined 
by a solution to (|j) and (|ic|), with 'P[+ lm = VlL lm = 0, for which Xi n (a n ) = 0. This approach 
is identical to that used by Kojima for the polar problem (Kojima 1993b). The idea seems 
consistent, but the definition of modes for a rotating star is a complicated issue that requires 
further discussion. 

The integration of (||) and (jjo|) is initiated at a small value of r using the regular power series 
solutions. In the case of © I use only the particular solution that corresponds to the coupling 
term Aj K m . The equations are then integrated to the surface of the star r = R. This part of the 
numerical calculation is straightforward. The main difficulty arises in the exterior of the star, 
and is associated with the well-known fact that the eigenfunction of a pulsation mode diverges 
as r —> oo. Specifically, a purely outgoing-wave solution will behave as Xi m ~ e lcrr * for large r. 
A pulsation mode that is damped with time according to an observer at fixed r corresponds to 
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Im a n < 0, and the desired solution to and © will therefore diverge as r increases. For a 
general complex frequency the decaying solution — that represents ingoing waves - will drown 
in the growing one and, unless Im er„ is relatively small, one cannot extract the asymptotic 
amplitudes x ou t and a,’i n from the numerical data. 

Several methods have been devised to deal with this problem (Kokkotas and Schutz 1992; 
Leins et al 1993; Andersson et al 1995). Most of these were taken over from the study 
of black-hole quasinormal modes. However, none of these methods is suitable for coupled 
equations. For example, while the WKB method of Kokkotas and Schutz (1992) can readily be 
used for a single second order ODE it is very difficult to use it for two coupled equations. In 
the slow-rotation problem one must deal with coupled equations. As long as the focus is on 
the spectral properties of the axial perturbations there is a way around this: By moving the 
term in the right hand side of © over to the left hand side and replacing X^ with X^ we 
only introduce an error of order e 2 . Thus, the new equation — that can be approached in the 
standard way — is consistent to O(e). This trick was used by Kojima in his study of the polar 
problem (Kojima 1993b), but I have chosen not to use it here. 

Instead, I have devised a “new” method to suppress the exponential divergence of the 
eigenfunctions. This method has the advantage that it will remain useful also for manifestly 
coupled equations. The integration of (|9j) and (|l^) is simply done using complex-valued r in 
the exterior spacetime. One would expect the exponential growth in X ~ e ± *°’ r ’* « e ±lar to be 
suppressed if the integration is performed along a straight line such that Arg(r) = —Arg(cr). 
Numerical tests show that this is indeed the case, and what remains is an (almost) purely 
oscillating function of r. This function can be integrated with satisfactory precision from the 
surface of the star out to a large value of |r|. The numerical solution is then matched to an 
expansion of form e ±lar * ^ a n (l/r) n at the final point of integration. This yields the values of 
the asymptotic amplitudes x ou t and x- m for any given complex frequency. Given these amplitudes 
one can iterate to find a frequency such that x- m {a n ) = 0, and identify a quasinormal mode of 
the star. 

To test the reliability of the new integration method I repeated the polar-mode calculations 
of Andersson, Kokkotas and Schutz (1995). Most of their mode-frequencies for a non-rotating 
polytrope could be reproduced with satisfactory precision. However, the new scheme is not able 
to get to the few ui-modes that are extremely rapidly damped. Moreover, the new integration 
approach is not accurate enough to enable iteration for very slowly damped modes, like the 
/-mode, in contrast to the method used by Andersson, Kokkotas and Schutz (1995). This is 
not a serious drawback, however. Slowly damped modes can be identified by a simple technique 
based on results for real frequencies [ the details of this method, which is an alternative to the 
standard resonance method (Chandrasekhar and Ferrari 1991b), will be described elsewhere]. 


4. The effect of rotation on the axial modes 

In the non-rotating case the symmetries of the perturbation equations dictate that the stellar 
pulsation modes be symmetric with respect to the imaginary a axis. That is, if a n corresponds 
to a mode then another mode is associated with — a n (the bar denotes complex conjugation). 
When the star is rotating this symmetry is broken. Each of the modes for a non-rotating star 
splits into 21 + 1 distinct modes, each corresponding to a different value of m (cf. the Zeeman 
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splitting in quantum problems), and the symmetry changes in such a way that <r nm is now 
associated with —a n _ m . 

The angular dependence of the perturbation functions is expressed in terms of spherical 
harmonics (cf. Kojima 1992). Specifically, the dependence on the azimuthal angle is e lmv . Since 
the time-dependence is assumed to be e~ lat positive frequency (Re a > 0) modes for m > 0 
will be prograde relative to the stellar fluid and m < 0 modes will be retrograde. Typically, 
retrograde modes will be destabilized by increased rotation, and the corresponding damping 
rates should decrease. For the prograde modes rotation will have a stabilizing effect, and their 
damping rate should therefore increase. The reason for this is the central mechanism of the 
Chandrasekhar-Friedman-Schutz instability (Chandrasekhar 1970; Friedman and Schutz 1978): 
Suppose that a mode has a frequency such that it is retrograde according to an observer in the 
fluid, but because of the dragging of inertial frames the mode appears prograde according to a 
distant observer. The distant observer — who measures energy according to his own time-line 

— finds that the mode carries positive angular momentum away from the star. But since the 
perturbed fluid actually rotates slower than it would in absence of the perturbation the angular 
momentum of the retrograde mode is negative. The emission of gravitational waves makes the 
angular momentum increasingly negative and leads to an instability. The same argument can 
be used to explain the destabilization of retrograde modes in general. 

A nonrotating perfect fluid star has no (non-trivial) axial pulsation modes in Newtonian 
theory. All the familiar fluid pulsation modes (the /, p and g modes) are associated with polar 
perturbations. But in the relativistic picture an infinite set of axial modes exist (Chandrasekhar 
and Ferrari 1991c; Kokkotas 1994). These are gravitational wave w-modes that exist because 
of the dynamic spacetime (see, for example, Kokkotas and Schutz 1992; Andersson et al 1996). 
The fundamental w-mode of a non-rotating neutron star typically has a pulsation frequency of 
the order of 10 kHz and damps out in a fraction of a millisecond. The modes become slower 
damped for more relativistic models. Furthermore, it is possible that an increased rotation rate 
will lead to relatively slowly damped (or even unstable) w-modes since retrograde modes should 
become slower damped as the rotation of the star increases. 

I have performed a sample of calculations for the polytropic equation of state 

p = 100km p 2 . (14) 

Figure [I] shows results for the slowest damped axial w-rnode of a stellar model with central 
density p c = 3 x 10 15 g/cm 3 . This leads to R = 8.861 km and M = 1.869 km = 1.27Mg, i.e. 
reasonable parameters for a neutron star. The results shown in Figure [j] nicely illustrate the 
effects of rotation on the stellar pulsation modes. As anticipated, each mode for the non-rotating 
star splits into 21 + 1 distinct ones. The prograde modes become shorter lived as the rotation 
rate increases, while the damping rate of the retrograde modes decreases. 

As already mentioned there are no non-trivial axial pulsation modes of a non-rotating star 
in Newtonian theory. But there is a set of Newtonian axial modes that are degenerate at zero 
frequency. As soon as the star rotates they become non-zero. These are the so-called r-modes 

— that are analogous to Rossby waves in the Earth’s oceans — and Newtonian calculations to 
first order in rotation show that they have frequencies of order f2 (Papaloizou and Pringle 1978; 
Saio 1982). Newtonian calculations also show that the r-modes remain degenerate at first order 
in fl. To unveil the detailed behaviour or the radial eigenfunctions and break the degeneracy 
one must proceed to higher orders (Smeyers 1980; Provost et al. 1981; Smeyers et al. 1981; Saio 
1982). 
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The r-modes have not previously been calculated in relativity, and the result of the present 
calculation is perhaps surprising: The relativistic analogue of the r-modes are unstable for all 
rotating stars. This result does, however, make perfect sense. A Newtonian calculation to first 
order in f l yields the r-mode frequency (Papaloizou and Pringle 1978; Saio 1982); 


oy = m£l 



KI+ 1 )) ’ 


(15) 


according to a distant observer. That is, these modes are retrograde relative to the fluid but 
will appear prograde at infinity. Such modes should (as already discussed) be destabilized by 
the emission of gravitational waves, and since the r-modes are initially (for a non-rotating star) 
situated at oy = 0 they should become unstable as soon as the star spins up. 

A sample of results for the r-modes and the polytropic equation of state 0 ) are shown 
in Table 0. The tabulated data are for the same stellar model that was used to generate the 
tr-mode data in Figure [3] (central density p c = 3 x 10 15 g/cm 3 ). It is interesting to note that the 
growth-times for these modes (represented by the “e-unfolding” time r e = 1/Im <j n ) are rather 
short. For example, for a millisecond pulsar the amplitude of the modes listed in Table |l] would 
multiply by a factor of three in a few seconds. 

The r-mode growth rates given in Table should, however, be viewed with some caution. 
The reason for this is the simple fact that, for small frequencies er ~ Q, the perturbation 
equations (9[) and (10) explicitly include terms that are formally of second order in e ~ Cl. As 
one can readily deduce from the data in Table [j] the numerical results scale with the rotation 
parameter as Im oy ~ e 2 as e —► 0. Hence, the estimated growth rates for the r-modes are due 
to the included second order terms — the eigenfrequency is real to first order in e. That is, the 
presented results for the imaginary parts of the eigenfrequencies are not obtained through a 
consistent calculation. It is possible (likely) that the results will change considerably once the 
calculation is extended to order e 2 . 

The numbers listed in Table [I] are nevertheless interesting. They provide the first suggestion 
that the r-modes are unstable in relativity. An argument based on the oscillation frequencies 
of the modes (which are of order e and should be obtained correctly in the calculation) and 
the mechanism behind the Chandrasekhar-Friedman-Schutz instability clearly establishes the 
existence of this instability. It is also far from clear how “wrong” the growth rates listed in 
Table |l] are. The error seems to be of the same order as the result itself, but this does not 
necessarily mean that the obtained result is useless as an order-of-magnitude estimate. What is 
clear at the moment is that a calculation that is consistent to second order in e, or an alternative 
estimate of the growth rate of the instability, is urgently needed. It will either show that 
Im oy ~ e 2 or establish that the true damping enters at higher orders. If the latter is the case a 
slow-rotation calculation of these growth rates seems very difficult. 

A second-order calculation would also be interesting since the degeneracy of the r-modes 
should be broken at order e 2 . Then r-modes with a different number of nodes in the radial 
eigenfunctions should lead to distinct eigenfrequencies (Saio 1982). 


5. Discussion 

I have presented the first numerical results for axial pulsation modes of a slowly rotating 
relativistic star. The calculation includes all terms of first order in e = R 3 /M. It includes 
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effects due to the coriolis force, but neglects effects due to the centrifugal flattening of the star. 
Such effects will appear at order e 2 and higher, and may well have a considerable effect on the 
various stellar pulsation modes. 

It has been shown that increased rotation has a stabilizing effect on prograde modes, while 
retrograde modes are destabilized. This was expected from the mechanism behind the familiar 
Chandrasekhar-Friedman-Schutz instability (Chandrasekhar 1970; Friedman and Schutz 1978). 
Specifically, it was demonstrated that rotation affects the axial gravitational wave w-modes in 
the anticipated way. 

I have also presented the first relativistic results for the so-called r-modes. The somewhat 
surprising result that these modes are unstable due to the emission of gravitational radiation 
even for the slowest rotating stars is interesting for many reasons. First of all, axial modes have 
not been considered in previous studies of stellar mode-instabilities (see, for example, Friedman 
1983; Lindblom 1986; Stergioulas and Friedman 1997). The present result shows that they are 
relevant, and that more detailed work is required in order for us to determine whether the 
unstable r-modes are astrophysically important. 

We need to establish if this new class of instabilities is stronger or weaker than the familiar 
instability associated with (say) the /-mode of a rapidly rotating star. The answer depends on 
the relative growth times of the various unstable modes. As already mentioned, the numerical 
results presented here (cf. Table [I]) should be used with considerable caution, but they will have 
to serve as a guideline until more reliable (higher order) data become available. 

If the unstable r-modes are to be of astrophysical interest they must grow on a time-scale 
much shorter than that associated with the damping due to viscosity, cf. the instability 
associated with the /-mode that is estimated to limit the rotation period of a 10 7 K neutron 
star at the level of 1-2 ms (Cutler and Lindblom 1987). To make a similar estimate based on the 
instability growth rates obtained in this paper is probably meaningless, especially since these 
results are expected to change (probably by several orders of magnitude) once more detailed 
studies are done. 

The r-mode instability is conceptually interesting, especially since it is active for ah, even 
the slowest, rotating stars. It is also possible that this new class of unstable modes of rotating 
relativistic stars is of considerable astrophysical relevance. The r-mode instability could, for 
example, limit the rotation period of nascent neutron stars. Whether this is the case or not 
is obviously unclear at present. It is a challenge for future work to provide an answer to such 
questions. 

It is a great pleasure to acknowledge enlightening discussions with (and/or helpful comments 
from) John Friedman, Yasufumi Kojima, Kostas Kokkotas, Norman Lebovitz, Lee Lindblom, 
and Bernard Schutz. I am also grateful to the Albert Einstein Institute in Potsdam for generous 
hospitality. 
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Table 1: r-modes for a T = 2 polytrope. This specific star has central density 3 x 10 15 g/cm 3 , which leads 
to R = 8.861 km and M = 1.869 km = 1.27Mq. All these modes are unstable due to the emission 
of gravitational waves. The dimensionless rotation parameter e, the rotation period 27 t/O, the pulsation 
frequency / and the “e-unfolding” time r e of the modes are listed. The data are for l = 2 and |m| = 2. 


e 

Period (ms) 

/ (Hz) 

Te (s) 

0.001 

404 

3.35 

1.36 x 10 5 

0.010 

40.4 

33.5 

1388 

0.025 

16.2 
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Fig. 1. - The effect of rotation on the slowest damped axial gravitational wave mode for a T = 2 
polytrope. This specific star has central density 3 x 10 15 g/cm , which leads to R = 8.861 km and 
M = 1.869 km = 1.27M©. The left panel shows the pulsation frequency of the star in kHz as a function of 
the dimensionless rotation parameter e = For this stellar model e « 0.2 would correspond to a 

pulsar with period 2 ms. The data corresponds to l = 2 and (from bottom to top) m = —2, —1, 0,1, 2. The 
right panel shows the e-folding time of the mode in ms as function of e. Here the curves correspond to (from 
top to bottom) in = —2, —1,0,1, 2. The m < 0 modes are retrograde while the m > 0 modes are prograde 
according to an observer in the stellar fluid. 




